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An H I shell-like structure associated with nova V458 Vulpeculae? 
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ABSTRACT 

We report the radio detection of a shell-like H I structure in proximity to, and probably associ- 
ated with, the nova V458 Vul. High spectral resolution observation with the Giant Metrewave 
Radio Telescope has made it possible to study the detailed kinematics of this broken and ex- 
panding shell. Unlike the diffuse Galactic H I emission, this is a single velocity component 
emission with significant clumping at ~ 0.5' scales. The observed narrow line width of ~ 5 
km s _1 suggests that the shell consists of mostly cold gas. Assuming a distance of 13 kpc 
to the system, as quoted in the literature, the estimated H I mass of the nebula is about 25 
Mq. However, there are some indications that the system is closer than 13 kpc. If there is a 
physical association of the H I structure and the nova system, the asymmetric morphology and 
the off-centred stellar system indicates past strong interaction of the mass loss in the asymp- 
totic giant branch phase with the surrounding interstellar medium. So far, this is the second 
example, after GK Per, of a large H I structure associated with a classical nova. 
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1 INTRODUCTION 

Classical novae (CNe) are interacting binaries, with orbital peri- 
ods of hours to around one day. A white dwarf (WD) accretes gas 
from a main sequence companion, and the build up leads to a ther- 
monuclear runaway (TNR) in the surface material. This generates 
a fireball, leading to the visual brightening that allows us to de- 
tect these events. This is accompanied by the ejection of around 
10~ 4 Mq of material at velocities of 200 to 5000 km s -1 (see e.g. 
IWarnej2008l for a recent review). These episodic stellar explosions 
are commonly detected at optical wavelengths and provide valuable 
understanding of the nuclear burning and accretion phenomena. 

Most of the CNe follow a maximum-magnitude-rate-of- 
decline (MMRD) relationship, where the time to decay by 2 or 3 
magnitudes from optical maximum (£2, £3 res pectively) is a reason- 
able i ndicator of the absolute magnitude (e.g. lDownes & Duerbeckl 
2000). As well as providing a distance estimate, this indicates the 
the underlying mechanism is fairly uniform across different exam- 
ples. The rate of decline also defines a speed class - from very 
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fast to slow - and it is evident that the fastest novae are also 
the most energetic (being intrinsically brightest and having the 
fastest ejecta). This suggests a link between the development of 
the ejecta, the progress of the explosion and the energetics (in- 
cluding any dependence on dredged up material) of the TNR (e.g. 
IShara. Prialnik & Shavivlll98rj|) . 

V458 Vul was identified as a CN on 2007 Augu st 8 (IAU Cir- 
cular 8861) with magnitude 9.5 dNakano et al]|2007l) . After reach- 
ing a maximum brightness of V = 8.1, the optical lightcurve 
showed unusual multiple peaks in the following weeks. If these ad- 
ditional peaks are ignored, it seems likely that for V458 Vul £3 ~ 
21 days, making this a fast nova. The MMRD relationship gives a 
distance of 13.5, 11.6 and 10.0 kpc to V458 Vul using t 2 , t 3 an d 
ti 5 respectively jDownes & Duerbeckl l2OO0l : IWesson et ai]|2008l) . 
iHenden & Munaril d2007l ) identified the star USNO-B1.0 1108 - 
0460444 as the "viable progenitor for the nova" base d on its colour 
and p roximity to the nova position (within 0".6; iNakano et al.1 
2007). Within 5" of the nova position, there is no other source in the 
USNO catalogue (completeness down to V = 21). Earlier spectra 
indicated V458 Vul to be an Fe II type nova, whilst later spectra 
were consistent with an He/ N type. This suggests t hat, in the spec- 
tral classifica tion scheme of IWilliams et all i 19941) . V458 Vul is a 
hybrid nova dPoggianil l2008). X-ray observations initially showed 
emission consistent with an expanding shock dTsuiimoto et al.1 
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Figure 1. Integrated H I column density map of the field showing large scale 
H I emission for the velocity range of Vlsr ~ "60 to —50 km s . The 
synthesized beam size is 30" X 30", and the colour scale is in mJy km s — 1 
per beam. The contour levels are for H I column density of (1, 2, 3, 4, 5) 
times 2.5 X 10 19 cm -2 . The lowest contour corresponds to a 3<r cutoff 
based on rms noise per channel (see [|2]for details). The optical position of 
V458 Vul is marked with a +. Note the broken shell-like structure north- 
east of the nova position. 



Table 1. Summary of the GMRT observation of V458 Vul 



Source position 


R.A. 19h54m24s.61 


(J2000) 


Dec. +20°52'52".6 


GMRT pointing 


R.A. I9h54m24s.3 


centre (J2000) 


Dec. +20°52'47".0 


Date of obs. 


June 11-12, 2009 


Bandwidth 


2.0 MHz, ~ 420 km s" 1 


No. of channels 


256 


Spectral resolution 


~ 1.6 km s 


Calibrators 


3C 48, 3C 286, 1924+334 


Observation time 


Total: ~ 10 hr, On source: ~ 4 hr 



2009), be fore a soft compon ent, reminiscent of RS Oph, came to 
dominate jDrake et alj|2008l) with a 12 hour oscillation becoming 
apparent, as this emission declined over 2008 (Swift nova consor- 
tium). In RS Oph this was associated with the unveiling of con - 
tinued hydrogen burning on the WD surface dOsborne et al.ll20llh . 
Follow up spectroscopy has established the the orbital p eriod of the 
binary system to be 98 min ( Rodriguez- Gil et alfcOld) . 

The field containing V458 Vul was observed both before and 
after the nova event , as part of the IPHAS Galactic Plane Survey. 
IWesson et al.l d2008l) identified a planetary nebula in the Ha im- 
ages, with a semi-major axis of 13". 5, with a position angle of 
~ 30 deg, centred to within 1" on the nova seen in the later IPHAS 
images. Further follow-up imaging and spectroscopic observations 
of the nebula showed an inner nebular knot, brightened by the flash 
ionization. This nebula, with an ionized mass of 0.2 M©, is sug- 
gested to be a 14000 years old planetary nebula orig inating from 
the c ommon envelope phase of the binary system dWesson et al.l 
120081) . 
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Figure 2. Zoom-in from FigurefT] showing the H I column density map of 
the broken shell. Contour levels and gray scales are the same as in FigurefT] 

In this paper, we report results of our radio observation at 1420 
MHz of the H I emission towards V458 Vul. This observation re- 
veals a shell-like H I structure close to, and probably associated 
with, the nova. The observation, analysis procedure and the re- 
sults are presented in Section [j2] Section |J3] contains discussion 
on possible implications of these results, and our conclusions are 
presented in Section ij4] 



2 OBSERVATION, ANALYSIS AND RESULTS 

V458 Vul was observed with the Giant Metrewave Radio Telescope 
(GMRT) with the aim of detecting any associated H I 21 cm emis- 
sion, and to search for possible radio continuum from the nova 
shell. Bandpass calibration for this observation was carried out 
with the flux and phase calibrators, but with frequency-switching 
at the first local oscillator, using a frequency throw of ±5 MHz. 
For details of GMRT o bservations using frequency switching, see 
iRov & Kanekarl J2007l) . Standard flux calibrators 3C 286 and 3C 48 
(at start and end of the observation, respectively) were observed for 
~ 30 minutes (10 minutes for each of the three frequency settings). 
A phase calibrator, 1924+334, was also observed in between the 
target source scans with each of the three frequency settings. A 
summary of the observations is given in Table[TJ 

All data were analysed in the Astronomical Image Process- 
ing System of the National Radio Astronomy Observatory (NRAO 
AIPS), using standard procedures. Some of the data were flagged to 
remove interference and instrumental misbehaviour. After remov- 
ing bad data, the standard calibration procedure was followed to ob- 
tain the antenna-based complex gains. The bandpass solution was 
derived from the ±5 MHz frequency setting scans of the calibra- 
tors, and the calibration was applied to the central frequency target 
source scans using interpolation between bandpass solutions. Next, 
the line-free channels were used to make a continuum image of the 
field. No radio continuum was detected at the nova position with a 
3cr limit of 1 mJy. 

The continuum subtracted residual data were imaged with 
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CLEANing in all channels. To recover weak and diffuse emis- 
sion structures, data from longer baselines were given a lower 
weight while making images. The restoring beam of the final im- 
age cube was 30" x 30" and the rms noise per channel was ^1.3 
mjy beam -1 . From this image cube, the moment maps were made 
using the channels with significant emission above 3<r noise level. 
Finally, this was converted to H I column density u nits assum- 
ing optically thin emission (Kulkarni&Heiles 1988). As shown 
in Figure [TJ we have detected discrete H I clouds in the vicinity 
of V458 Vul with velocities ~ —55 km s~ . If the gas velocity is 
due to Galactic rotation only, the kinematic distance of the struc- 
ture is 13.5 — 14.5 kpc. A broken clumpy shell-like H I structure 
north-east of the nova position is clearly visible in this image. No 
such dumpiness at similar scales was seen in the H I maps at the 
positive velocities, indicating that the H I emission component is 
smoothly distributed over large angular scales, which is typical of 
the hot atomic gas detected in the H I 21 cm emission line in our 
Galaxy. The details of column density and velocity field of the bro- 
ken H I shell are shown in Figure [2] and Figure 0a) respectively. 
Figure [3jb) shows the H I emission in different channels for the 
velocity range of —46.8 to —60.0 km s _1 , covering the spectral re- 
gion of interest. The central channels show a line of sight velocity 
difference of ~ 3.5 km s _1 between the south-east and north-west 
rim. 

The integrated spectrum of the H I shell is shown in Fig- 
ure |3]c). The blue line is the best fit Gaussian profile with a full 
width at half maximum (FWHM) of only 4.9 ± 0.3 km s" 1 . If 
there was no non-thermal broadening, this would correspond to a 
temperature ~ 525 K. However, as noted above, at least ~ 3.5 
km s _1 is due to the internal dynamics of the H I structure. There 
may also be some contribution from turbulence. So the actual ki- 
netic temperature is likely to be considerably lower. We also used 
the line profile to estimate an H I mass of (25.0 ± 2.3) x df 3 Mg, 
where di3 is the distance in units of 13 kpc. 

The H I emission towards V458 Vul (in the first Galactic quad- 
rant) detected near Vlsr = ^53 km s" 1 is a single narrow emis- 
sion component while rest of the H I towards V458 Vul (and also 
towards the nearby phase calibrator 1924+334) at positive veloc- 
ities has multiple components spread over > 50 km s _1 . The H I 
emission near —53 km s _1 shows clumping at ~ 0.5' scales. In 
contrast, the H I emission detected at positive velocities is more 
widespread and of relatively uniform brightness over the GMRT 



primary beam with no significant clumping below ~ 10'. The nar- 
row H I line width of the shell suggests that the gas is cold. How- 
ever, the estimated size and mass (~ 7.6 x di3 pc and 25 x d\ 3 
Mq respectively) are larger than the typical size and mass scales 
(1 — 2 pc and 2 — 5 M0 respectively) of the Galactic cold neutral 
medium. These characteristics indicate a possible peculiar origin 
(e.g. association with V458 Vul, shaped by an early phase of mass 
loss and interaction) of this shell. 



3 DISCUSSION 

In this paper we report the discovery of a large asymmetric shell- 
like structure in H I to the north-east of V458 Vul. At the maxi- 
mum likely distance of 13 kpc, the average projected linear diame- 
ter of the H I structure is about 7.6 pc. For comparison, GK Per 
has a planetary nebula of size 6 pc with associated H I struc- 
ture dSeaquist et al.ll9 89: Anupa ma & Kant haria 2005). Moreover, 
there have been reports of planetary nebulae observed in dust emis- 
sion at infrared wavelen gths which are several parsecs in size 
dWeinberger & Arvali r2004). So, a size of a few parsecs is not in- 
consistent with our interpretation of this structure as being asso- 
ciated with V458 Vul. V458 Vul is thus the second classical nova 
found to be associated with a large H I structure. 

If the H I structure and the nova have a physical association, 
we may interpret this shell as a remnant of the interaction with 
the surrounding interstellar medium (ISM) in the asymptotic gi- 
ant branch (AGB) phase of evolution of the white dwarf. The H I 
structure, which appears to be ram pressure swept gas, has a size 
of about 2'. 4 x l'.8 at a position angle of ~ 30deg and a sys- 
temic velocity of —53.2 km s _1 . Interestingly, the position angle 
is very close t o the orientation of the planetary nebula detected by 
IWesson et all fc008t) . though the angular size of the planetary neb- 
ula is much smaller (~ 27"). This, we believe, supports a common 
origin to this axis. It might indicate the ISM density structure or the 
axis of mass loss in the binary system. 

3.1 The asymmetric H I shell: evolution in the AGB phase 

The asymmetric nature of the structure and the off-centred stel- 
lar system may be attributed to a number of causes includ- 
ing the effect of the interaction of the stellar wind with the 
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ISM, the eff ect of the binary stellar sy s tem and magnetic field 



effects (e.g. iKwokl Il982l: ISokerl Il98s j : Ichevalier & Luol 



Dwarkadas. Chevalier & Blondinl 1 19961 : iBIackman et al 
Soker & Rappaportl feoOlh 
Simulations 



1994 



2001 



Borkowski. Sarazin & Sokerl dl99d) , 
Soker. Borkowski & SarazuT dl99lh and more recently by 
Wareing. Ziilstra & O'Brienl J2007I) have examined the interaction 
of ambient ISM with systems with continuous or episodic mass 
loss. Though these are mainly in the context of planetary nebulae, 
they can explai n several of the observed features in this asymmetric 
H I structure. IWareing. Ziilstra & O'Brienl d2007h . in particular, 
show that, in about 3 x 10 4 years into the post-AGB evolution, the 
central stellar system is not at the centre of planetary nebula and 
the nebula is highly asymmetric with the brighter emission in the 
dow nstream regions. For a density of 2 cm -3 of the surrounding 
ISM. IWareing. Ziilstra & O'BrienM2007t) find that atomic densities 
up to 50 — 100 cm~ s can be present in the ram pressure stripped 
flow in the downstream region. This is similar to and consistent 
with our estimated density of about 10 cmT A in the H I structure 
for a distance of 13 kpc. 

Severa l AGB stars have been observ ed and detected in H I 
emission dGerard & Le Bett re 2006, 2007) and many of these stars 
show a circumstellar envelope of size ~ 1 pc. More recently 
iMatthews et al.l d201 ll) report the detection of a circumstellar wake 
and a compact high ve locity cloud near the AGB star X Herculis. 
iMatthewset al. (2008) reported the detection of an H I tail coin- 
cident with part of t he 4 pc long t ail in the far ultraviolet trailing 
the AGB star o Ceti. lMartin et all d2007h suggest that the long tail 
in this case is a result of the large space velocity (~ 130 km s _1 ) 
of the Mira variable (which is a binar y system) and the interaction 
betwe en the ISM and the stellar wind. War eing. Ziilstra & O'Brienl 
d2007l) concluded from their simulations of the system, that the long 
tail is a result of mass loss from the AGB star over 4.5 x 10 yr. 

Thus, we infer that (i) large H I structures and tails as- 
sociated with stellar systems in the AGB or post-AGB phase 
are fairly common, and (ii) the observed asymmetry of the H I 
structure associated with V458 Vul may be due to the motion 
of the nova system with respec t to the ambient ISM (see also 
iBode. O'Brien & Simpso"r]|2004 for evidence of this in GK Per). 
Unfortunately, such structures are difficult to detect due to the pres- 
ence of Galactic emission kinematically corrupting the circumstel- 
lar emission and the inherently weak signal which limits detections 
only to nearby systems with peculiar velocities. 



3.2 Distance, age, and kinematics 

Distance estimates to novae have always been difficult to obtain 
and V458 Vul is no exception. From the different estimates in liter- 
ature, it is clear that the dist ance estimate t o V45 8 Vul is uncertain 
by at least a factor of two. Wesson et al. ( 2008) estimate a range 
of 10 to 13.5 kpc from MMRD relation. Their estimations of the 
distance, based both on the light travel time argument and on the 
assumption that mean radial velocity of Vlsr = —60.6 km s" 1 
is due to Galactic rotation , are also consistent with a distance of 
~ 13 kpc. Poggiani (2008) estimate the distance to be in the range 
6.7 to 10.3 kpc usi ng several methods, some o f them similar to 
IWesson et all j2008). Recently, Rajabi et al. l i2012l ) have reported a 
distance of 9.9 — 11.4 kpc by modelling the expansion of ejecta 
based on optic al interferometric observation. On the other hand, 
iGordonl f201lh estimates the distance to be 5 to 6.5 kpc using the 
metho d of light echoes (although light echo modelling bv lHounsell 
I2012L gives d ~ 13 kpc). Note that the line of sight velocity may be 



Table 2. Derived physical parameters of the H I structure for different as- 
sumed distances. The columns in the table are (1) Distance to the nova (D), 
(2) Size of the H I shell (dni), (3) H I mass, (4) Estimated dynamical age 
(^dynamical)* (5) proper motion (p.) corresponding to the observed separa- 
tion of V458 Vul from the centre of the H I shell, and (6) Tangential velocity 
(vt) corresponding to a proper motion of 20 mas yr~ 1 . See i]3.2l for details. 



D 


dm 


M H i 


^dynamical 






kpc 


pc 


M 


years 


mas yr 


km s 


13.0 


7.6 


25.0 


3.7 X 10 5 


0.27 


1200 


6.5 


3.8 


6.3 


1.8 x 10 5 


0.56 


600 


5.0 


3.0 


3.7 


1.4 x 10 5 


0.71 


474 


1.0 


0.6 


0.15 


2.8 x 10 4 


3.57 


95 



anomalous (like the high velocity clouds), and the actual distance 
may be very different from the kinematic distance. Also, with no re- 
ported observation before August 8, both the maximum magnitude 
and the rate of decline is highly uncertain in this case. Moreover, 
multiple peaks in the declining part of the optical lightcurve makes 
defining ti and tz, and hence using MMRD relationship, non-trivial 
in this case. It is also important to note here that, even with ex- 
act determination of peak apparent magnitude and rate of decline, 
the intrinsic scatter of the MMRD relation (e.g. ldella Valle & Livid 
1995; Downes & Duerbeck 2000) and the uncertainty of extinction 
can result in a factor of ~ 2 uncertainty in the derived distance. 

Obtaining a reliable distance estimate to the nova system is 
important in determining the physical parameters of the H I struc- 
ture. In Table© the derived parameters for the system are presented 
assuming different distances to the nova. The columns in Table [2] 
are (i) assumed distance, (ii) size and (iii) mass of the H I shell, (iv) 
dynamical age of the shell for an average radius of 60" assuming 
an expansion velocity of 10 km s" 1 , (v) proper motion correspond- 
ing to the observed separation of V458 Vul from the centre of the 
H I shell, and (vi) tangential velocity corresponding to a proper 
motion of 20 mas yr -1 (see more on this below). Note that even if 
the line width is ~ 5 km s _1 , as we do not know the inclination, 
we have used a typical value of 10 km s~ 4 for expansion veloc- 
ity. If the expansion velocity on the sky plane is the same as the 
line of sight expansion velocity, the age will be about factor of two 
higher. Interestingly, in spite of the significant spatial displacement 
between the H I structure and the planetary nebula, and the differ- 
ence in their estimated ages, their radial velocities are reasonably 
close (Vlsr = —53.2 and —60.6 km s~ 4 respectively). 

The star USNO-B1.0 1108 - 0460444 has be e n iden tified 
as the progenitor of V458 Vul by Irlenden & Munaril d2007n . The 
proper motion of the corresponding star listed by NOMAE0 is 
about — 12±5 mas yr -1 in Right Ascension and 16±2 mas yr -1 in 
Declination (with the highest value of the probability estimator as- 
signed for the likelihood that the quoted proper motion is correct). 
This is roughly 20 ± 3.5 mas yr~ in the north-west direction. For 
a distance of 13 kpc, the proper motion translates to an unreason- 
ably large value of tangential velocity of 260 ± 44 AU per year i.e. 
1200 ± 200 km s" 1 (see Table©. On the other hand, the required 
velocity, to displace the nova by the observed separation of ~ 100" 
from the centre of the H I shell over this dynamical timescale, is 
only ~ 17 km s _1 or 0.27 mas yr -1 at 13 kpc (see Table©- This 
is the vector sum of the proper motion of the star and the ejecta 
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(which may be misaligned; Ge rard. Le Bettre & Libertll201 lh . So, 
the proper motion of V458 Vul is expected to be smaller than this. 
Also, V458 Vul is centr ed to within l" of the 14000 years old plan- 
etary nebula dWesson e t al. 2008), implying proper motion smaller 
than 0.07 mas yr _ , It is more likely that the NOMAD proper mo- 
tion is spurious for such a faint object (e.g., due to issues with the 
older plates used to derive proper motion). 

We also examined the possibility that USNO-B1.0 1108 - 
0460444 is not the progenitor of V458 Vul. At a distance of ~ 13 
kpc, such a progenitor must be fainter than V = 21 (My > 3.5 
for Ay = 1.95). Taking a conservative lower li mit of My = 10.5, 
and using a typical stellar luminosity function dKroup a 2002), the 
average number density of stars with 3.5 < My < 10.5 is 
~ 2.5 x 10~ 2 pc -3 . Considering the position uncertainty to be 
~ 0".6, the volume, defined by a circle of 0".6 radius (on the 
plane of the sky centred at the nova position) and a radial distance 
of 13 ± 2 kpc (2 kpc centred at a distance of 13 kpc), is about 
9 pc 3 . Within this volume, the average expected number of pro- 
genitor stars with the right My is about 0.22 (0.62 for 1" radius). 
The associated Poisson distribution probabilities are 0.18 and 0.33 
for 0".6 and 1" respectively. Clearly, the possibility of an alterna- 
tive nova progenitor at a distance of ~ 13 kpc is small. 

If the H I structure has a physical association with the nova, 
as implied by various pieces of evidence above, then an indirect 
argument in favour of a smaller distance to the nova is from the es- 
timated mass of the H I shell. At a distance of 13 kpc, the enclosed 
H I mass is ~ 25 Mq which is too large to be explained as swept- 
up mass. If the distance is < 6.5 kpc, the H I mass may originate 
from a massive super- AGB progenitor. Such pro genitors evol ve to 
either neutron stars or massive ONeMg WDs dHerwig|2005l) . In- 
terestingly, the massive (> 1 Mq) WD in V458 Vul is consistent 
with this scenario. If the NOMAD proper motion, as argued above, 
is incorrect, based on parameter values in Table[2] a distance < 6.5 
kpc seems plausible. It is important to constrain the proper motion 
of the faint nova progenitor, and if possible find tracers from the 
nova outburst which would give independent distance estimates, to 
cross check this suggestion of a smaller distance to V458 Vul. 



3.3 Possible large scale H i structure 

The above discussion has been mainly restricted to the shell-like 
expanding structure northeast of V458 Vul. However high veloc- 
ity H I has been seen over a much larger region around V458 Vul. 
There are discrete clouds scattered around the nova with the central 
few arcminutes being devoid of H I (see FigureQj. One possibility 
is that there is another larger shell which traces the mass loss his- 
tory at a still earlier evolutionary epoch of the nova binary system. 
This large shell (approximately 8'), with longer axis almost parallel 
to the smaller H I shell we discussed above, also has the nova off- 
centred. Moreover, this large shell seems to show weak evidence 
of expansion with the north/north-eastem parts moving at —49 
km s _1 and the southern parts being close to —55 km s . How- 
ever this needs confirmatory observations especially since GMRT 
is expected to be missing flux on larger angular scales. No clear 
correlation is seen between the H I and the IRAS 60 and 100 /J,m 
dust emission, though the H I shell appears to be located close to 
the edge of a dust cloud. 



4 CONCLUDING REMARKS 

From the above discussion, we conclude that we have detected an 
expanding shell of cold H I possibly associated with V458 Vul. If 
this is the case, it appears to be the remnant of the interaction of the 
nova system with the ISM in its AGB phase of evolution, as indi- 
cated by its slowing down and its peculiar morphology. However, 
at a distance of 13 kpc, the H I structure would have an unreason- 
ably high mass of 25 Mq . There are some indications that the nova 
is closer than 13 kpc, and further independent constraints on the 
distance to the nova will be very useful. 

Finally, it is relevant to mention that a large number of such 
systems might be distributed in the Galaxy. While we have detected 
the H I shell associated with V458 Vul, it is likely that there are 
several such structures which have not been detected in H I emis- 
sion either because of the inadequate spectral resolution of such 
observations, or because they are so far from the central stellar sys- 
tem that the association is not obvious. More observations targeting 
such systems may be useful to understand the interaction with the 
surrounding ISM and the mass loss phase of the evolution of clas- 
sical novae progenitor system. 
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